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Abstract 

In this work we find cosmological solutions in the brane-bulk system starting from a 5-D line 
element which is a simple extension, for cosmological applications, of the pioneering Randall- 
Sundrum line element. From the knowledge of the bulk metric, assumed to have the form of plane 
waves propagating in the fifth dimension, we solve the corresponding 4-D Einstein equations on 
the brane with a well defined energy-momentum tensor. 
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1. Introduction 



Brane cosmologies are often studied in the framework of five-dimensional (5-D) Einstein 
equations in the bulk. The effective 4-D gravitational equations in the brane without cur- 
vature correction terms were first obtained by Shiromizu, Maeda and Sasaki [1]. However, 
even taking more generalized gravitational actions, the derived 4-D Einstein equations do 
not in general form a closed system due to the presence of a Weyl term which can only be 
specified in terms of the bulk metric, so other equations are to be written down and different 
procedures arise in splitting the non-Einsteinian terms between bulk and brane [2] . Because 
the specific form of a solution in the bulk is in general rather cumbersome due to the number 
of terms and parameters in the equations, we shall consider, as a guide for further work, a 
model simple enough to allow obtaining not trivial exact solutions. This paper is organised 
as follows. In the next Section we summarize the effective brane equations obtained by Ko- 
finas [3] when the intrinsic curvature scalar is included in the brane action. In Section 3 
we transform the static line element proposed by Randall and Sundrum [4] into a dynamical 
one containing a three-space with constant curvature and obtain, using the 5-D Einstein 
equations, the corresponding energy-momentum tensor in the bulk. In Section 4 we find 
the related bulk metric in coordinate systems commonly used in cosmological applications. 
Solutions in the brane, assumed infinitely thin and Z 2 symmetric in the bulk, are found in 
Section 5. Finally, conclusions are given in Section 6. 

2. Braneworld Einstein field equations 

In this section we recall the effective brane equations obtained by Kofinas [3], which we 
shall use in the following, giving a brief account of their derivation. Once we have solved the 
equations in the bulk, the form of the induced equations will allow us finding brane solutions 
following the methods of General Relativity with a well defined energy-momentum tensor. 
The starting point in [3] is a three-dimensional brane E embedded in a five- dimensional 
spacetime M. For convenience a coordinate y is chosen such that the hypersurface y = 



2 



coincides with the brane. The total action for the system is taken to be 

+ / V-^gLr* d 5 x+ [ ^/-WgL™ at d 4 x (1) 

J M JT, 

Bulk indices will be denoted by capital Latin letters and brane indices by lower Greek letters. 
Varying (1) with respect to the bulk metric qab one obtains the equations 

< 5 >G* = -A 5 5 B + 4 (< 5 >Tf + (ioc) T*8{y)) (2) 

where 

- 2 J^^G B A -ni^n +^h B A ) (3) 

is the localized energy-momentum tensor of the brane. ^Gab and ^Gab denote the Ein- 
stein tensors constructed from the bulk and the brane metrics respectively, while the tensor 
Hab = Qab — nAfiB is the induced metric on the hypersurfaces y = constant, with n A the 
normal unit vector on these 

rA 

n A = -^, n A = (0,0,0,0,$) (4) 

The way the coordinate y has been chosen allows to write the five- dimensional line element, 
at least in the neighborood of the brane, as 



(loc)rpB _ _j_ 
J- A " 



dS 2 = g AB dx A dx B = 9ilv dx»dx v + ^dy 1 (5) 

Using the methods of canonical analysis [5] the Einstein eqs. (2) in the bulk are split into 
the following sets of equations 

K»K; -K 2 + (4) i? = 2 (A 5 - 4 (5) 7*) (6b) 
+ $ K K» - $ + g» x $ ;Al/ = - k\ $ (t° c )T» - 1 (loc ^T 5{y) - k\ $ (5) T^ 

- 4 $ (5) T^ + \$(k 2 5 (5) T - 2A 5 ) 6% (6c) 



dy 



where K„ v is the extrinsic curvature of the hypersurfaces y = constant: 



The Israel conditions [6] for the singular part in eqs. (6c) are 

[K*\ = -k 2 5 $ (Ve) T M _ 1 (ioc) T ^ 

where the bracket means discontinuity of the quantity across y = and $o — *&(y — 0). 
Hereafter, considering a Z 2 symmetry on reflection around the brane, eqs. (3) becomes 

= -A 4 # + + — (K1 - K8») (9) 

where = K^{y = + ) = — Kjf(y = 0~) and r c = k\/ 'k\ is a distance scale. Equations (9) 
assume the form of usual Einstein equations with the energy-momentum tensor splitted into 
the common brane energy-momentum tensor plus additional terms which are all multiplied 
by l/r c . The tensor satisfies the usual conservation law ^Tjf = provided = 0, 
which means no exchange of energy between brane and bulk. The quantities are still 
undetermined, however additional information can be obtained from the geometrical identity 

^R A BCD = ®R% KL hi h N B h K c h L D + (K$ K BD - K% K bc ) (10) 

Taking suitable contractions from the above relation it is possible to construct the four and 
five-dimensional Einstein tensors and, making use of the bulk Einstein equations, to get 
finally the parallel to the brane equations 

+ (kX - K" x Kt) + \ (KlK X K - K>) 8» - sT ^C y Kyu (11) 

Here is the "electric" part of the bulk Weyl tensor, while T and C are the limiting 

values of those quantities at y = + or CT. Now, equating the right-hand sides of the inde- 
pendent eqs. (9) and (11), one gets an algebraic equation for which can be substituted 
in (9) once the equation has been solved. However the system of Einstein equations for the 
brane metric so obtained is not, in general, closed due to the presence of the Weyl term, so 
one has to solve the Einstein equations in the bulk in order to determine the Weyl tensor on 
the brane. This one will be the method followed in this paper. A different approach, which 
instead does not assume a bulk geometry, starts from deducing a brane dynamics and then 
searches for a bulk geometry in which the brane can consist its boundary [7]. 
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3. The model 



Following a possibility suggested in [8] we shall consider 5-D Einstein equations strictly 
in the bulk, i.e. without the brane energy-momentum tensor with its delta distribution, so 
when searching, later on, for solutions in the brane we shall have to take limiting values of 
bulk quantities as discussed in the context of eqs. (11). As a simplifying device for dealing 
with these equations, we select a five- dimensional line element with reasonable physical 
properties, calculate the corresponding energy-momentum tensor and, once solved 

eqs. (11) for the brane metric, we can obtain, from eqs. (9), the 4-D cosmological constant 
A 4 and the effective energy- momentum tensor on the brane. More in detail, let us consider 
the static Randall-Sundrum line element 



d S* = e -2*v [A 2 da 2 k -dt 2 )+dy 2 



(12) 



containing a three-space with constant curvature. Here k and A are, respectively, the con- 
stant scale factors for the extra dimension y and for the ordinary three-space and da\ is the 
line element on maximally symmetric three-spaces with curvature index k = +1, 0,-1: 

dr 2 



da\ 



+ r 2 (<M + sin 2 tfcV) 



(13) 



1 — kr 2 

Since our purpose is to describe the time evolution on the braneworld, we need to transform 
the static bulk solution (12) into a dynamical one. We follow the procedure used in [9,10,11] 
where dynamical solutions are derived from the static Randall-Sundrum metric by gener- 
alized boosts along the fifth dimension. Applied to the actual case, we write the required 
transformations as 



1 



,-xt 



t 



X 



X_ 

K 2 



X 



(14) 



Ky 



e Ky = 



+ e 



-xt 



K z 



where k = 



K 



X = 



X 



X 



K 



and x is a constant responsible of the boost in the 



1 - 



X 



K 



(t, y) spacetime. 



As a result the metric (12) becomes, dropping the bar: 



L_-(i 2 ^-e- 2 ^^ + e 2 -^) (15) 



dS z = 7 _ 
with 1 2 = (1-^)A 2 . 

Having specified the metric components and consequently the components of the Einstein 
tensor, we have 

G; = G% = G* = -6(?-J?)-- s £-- (16a) 

GHG; ._ 6(? _, V _JA_ (16b) 



where 



a (£ y ) = = (17) 

v ' i " e -xt _|_ e ny _ 1 v > 



By comparison with (16) and taking into account that G r t and G\ are both equal to zero, 
we obtain 

A 5 = 6 (£2 _ ~2 } = _ g ^ (5) T B = diag(p; p> p> _ p? p±) (lg) 

with 

2 ^ 2 2 3 /c , . 

a 2 (t,y) a 2 (t,y) 
It is worth noticing that the bulk sources contain, besides the cosmological term which can 
be interpreted as proportional to the pressure and to the density of vacuum fluctuations 
[12], terms proportional to the curvature index k which call in mind, even if k and a 2 (t,y) 
refer to different spaces, a relationship with the "pressure" and "density" coming from the 
curvature of the universe [13]. The line element (15) is a particular realization of the metric 

dS 2 = a 2 (t, y) da\ - n 2 (t, y) dt 2 + $ 2 (t, y) dy 2 (20) 

commonly used in cosmological applications, so it is worth obtaining the above metric coef- 
ficients solving Einstein with a generic cosmological constant A 5 and the tensor whose 
components are given in (19). 
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4. Solutions in the bulk 



In the coordinate system (20) the non-vanishing components of the Einstein tensor G\ 



are 



Gr (~fQ /-up 
r — W — U <p 



$ 2a $ (2a h 
$ a $ 1 a n 



a ( a 2h\ 
a la n J 



G y 

y 



G: 




(21a) 
(21b) 
(21c) 
(21d) 



Here a dot and a prime denote partial derivatives with respect to t and y, respectively. 
The specific form of the solution in the bulk is in general rather cumbersome, but some 
simplifications can arise if one makes the assumption of plane waves solutions, i.e. if one 
assumes that the metric coefficients are functions of the argument u = (t — Xy) or of 
v — (t + Xy) [14,15,16,17]. In the following we concentrate on solutions which depend on 
(t-Xy): 

a = a(t-Xy), n = n(t-Xy), $ = $(£- Ay) (22) 



From the equation G* = we get 



constant 



(23) 



where the superscribed asterisk * denotes derivative with respect to u. Now considering the 
history of the three-brane as described by a point trajectory in the (t, y) spacetime it is 
possible to introduce a Gaussian normal coordinate system where $ = 1. Another ansatz 
where n — 1 is made in the literature [18,19] so we shall separately consider these two 
additional assumptions to solve the Einstein equations in the bulk. 
Let us begin with the choice n — 1. One has from (23) 



. a 
$ = - 



(24) 
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with P a constant. Now G\ = G v y and the Einstein equations reduce to 



2 ' 

— h4oo + S = A 5 a 2 + I3 2 \ 2 (25a) 

a 

a * a * + a 2 =^a 2 + p 2 \ 2 (25b) 

Subtracting (25b) from (25a) we obtain an equation which is the first derivative of (25b) 
with respect to u, so both eqs. (25) are satisfied solving 

(a 2 ) = -A 5 a 2 + 2p 2 \ 2 (26) 
3 

The solution is 



3 2 \ 2 

a 2 = c 1 sinhy^A 5 (£- Xy) + c 2 cosh ^ | A 5 (t - Xy) - 3 ^— (27) 

with ci and C2 suitable constants. Requiring a 2 (0) = and putting for future simplifications 
ci = we have 

a 2 = ^sinh 2 y5 (t _ Ay) (28) 
The other metric coefficient is 



^ = A 2 oo8h^^(*-Ay) (29) 

* 

We notice that the choice ci = implies (a 2 )(0) = and leaves the possibility of taking 
either sign for the 5-D cosmological constant. 
When $ = 1 one has from (23) 

* 

n = - (30) 

a 

with a a constant. One has again G\ = G y y and the Einstein equations reduce to 

r\ * =K * . 

a a ** *2 A 5 2 cr 

— j \-4aa+a = - — a +— (31a) 

a A A 

A 2 

** , * 2 A 5 2 , a / Q1 n 

flfl+fl =- TT2 a + ~vl (Sib) 
o A A 

Subtracting (31b) from (31a) we obtain an equation which is the first derivative of (31b) 
with respect to u, so both eqs. (31) are satisfied solving 

(5) = -^ + 2?! (32) 
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The solution is 

1 / 2 . , i 1 / 2 . , . x 3 a 2 



a 2 = ci sinh^y- | A 5 (i - Xy) + c 2 cosh-^- | A 5 (t - Xy) - (33) 
Requiring a 2 (0) = and putting for future simplifications C\ = we have 



A 5 A V 6 

The other metric coefficient is 



a 2 = -^sinh 2 ij-^(t-Ay) (34) 



n 2 = loosh 2 i^-^(t-Ay) (35) 

* 

Here again the choice ci = implies (a 2 )(0) = and leaves the possibility of taking either 
sign for the 5-D cosmological constant. When Z 2 symmetry on reflection around the brane 
at y = is assumed, the extra-coordinate y in the above solutions must be replaced by \y\. 

5. Solutions in the brane 

The 4-D line element can be written as 

ds 2 = al(t) da\ - n\{t) dt 2 (36) 

where the subscript o here does not mean "calculated putting y = in the corresponding 
bulk quantities" because those metric coefficients were found without taking into account 
the matter content of the brane, but means "calculated solving Einstein equations (11) on 
the brane with (^G^ derived from the line element (36)". Before searching for plane wave 
solutions in the brane let us return, for completeness, to the case treated in Section 3 with 
5-D line element given in (15). The relevant Einstein equations obtained from (11), are 

(4W _ (4)™? _ (4W - _ k °o , 2 «o^o _ 2 "o _ _ A 5 _ 7 ~ 2 _ ke~ 2xt 

a «o n o a o n o a o n o z A z 

3 hp" 2 * 1 

(37b) 



dg ! 2d n 
agn 2 , a ng 


2 d 


A 5 
2 


7£ 2 - 


a 2 


3d 2 


A 5 
2 


7£ 2 - 



A 2 



where A 5 is given in (18). 
The solutions are: 



a 2 = Cl e 2 ^=A 2 e 2 ^ (38) 

n 2 Q = - 2 ^0 = 1 (39) 

|(A 4 + K 2 p)- fc 



where A 4 and p, written explicitly below, are respectively the cosmological constant and the 
matter density in the brane. With the choice c\ = A 2 the brane metric can be obtained 
evaluating the corresponding bulk metric at y — 0. Comparing eqs. (37) and (9) we have 
defined the 4-D cosmological constant, including a possible contribution coming from the 
brane tension, as 



A 4 = ^ + 7£ 2 = (40) 



A 5 _4^+3/ 

~o~ t ' K ~ 2 

K 

The remaining terms are the components of the effective energy-momentum tensor of the 
brane. This can be considered as the stress tensor of a perfect fluid at rest, because G rt = 
in the frame given by (36), with pressure p and density p given by 



3A; 

3 /t 2 Oq 



P = ~oP, P = ZTZ2 ( 41 ) 



Here and in the following we shall attribute to a "fluid" also quantities proportional to the 
curvature index k and to the cosmological constant A 4 ; moreover all the fluids taken into 
account may cause violations of some of the energy conditions. Comparing the values of the 
Hubble constant obtained from its definition and from the Friedmann equation (37b) 



a n J 3 % 3 3 

we obtain x — H — a/A 4 /3 for any value of the curvature index k so a oc e Ht as in the de 
Sitter model, being now t the proper time in the brane. Clearly the deceleration parameter 
q = — (a d )/a 2 and the density parameter = A 4 /(3iJ 2 ) have the values q — — 1 and 

n A = i. 

The next case we shall treat corresponds to the choice n — 1 in the bulk. The relevant 
Einstein equations now are 

k Oq 2 d n 2 a A 5 2 A 5 k 



a 2 a 2 n 2 a n 3 a n 2 2 /XT 6/? 2 A 2 . , 2 /A 5 

3 smh W — t — ; smh \ — t 

V 6 A 5 V 6 

3k 3al A 5 A 5 3 k 



V 6 A, V 6 



(43b) 
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The solutions to the system (43) are 

2-k/(k+/3 2 X 



ai= Cl [sinh^/^) (44) 

nl= - 2 % (45) 

-| (A 4 + k\ (pi + P2 )) - k 

Comparing eqs. (43) and (9), we have defined the 4-D cosmological constant as A 4 = ^ 
and have considered the remaining terms as the superposition of two different fluids with 
pressures and densities given by 

2A4 -2(fc+/3 2 A 2 )/(fc+2/3 2 A 2 ) 



Pi — - g Pi) Pi — t^=- a I 4 " J 

k\ sinh 2 J^t 

p, = -\p,, P- 2 = „ 2 , * ^ oc ^-^^ (47) 



A 4 V 3 

From (43b) we obtain the Friedmann equation 

\a Q noJ 3 3 

Here the density p = p\ + P2 appears linearly as in the usual Friedmann equation. This 
linear dependence is common to other induced gravity models in the literature [20] and 
appears, of course with different values of the quantities involved, in all the brane solutions 
found in this paper. Also we notice that in the actual case it is neither possible to recover 
the brane metric evaluating the bulk metric at y = nor in general to give the explicit 
dependence of the scale factor on the proper time r. As a qualitative estimate one can say 
that Oq, starting from zero, either increases indefinitely or oscillates depending on the values 
of the parameters which appear in eq. (44). Exact results can however be obtained if we 
assume k — 0, which is a value strongly suggested for our universe at the present time. For 
a spatially flat universe we have from eq. (48) 

a oc sinh 2 A /5 T; #2 = A4 C0th 2 /A7 r d \ ( x M (4g) 

V 12 ' 3 V 12 y 2 V 3#V V ; 

Recalling that the density parametr fi A = A 4 /(3iJ 2 ) varies in the interval [0,1], one ob- 
tains for the deceleration parameter — 1 < q < — \ and for the age of the universe 
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t = 2/H ^a/T7(^a)o cotfT 1 ^/V(^a)o) where the subscript here refers to present values. 
While the range of values for q is in agreement with the observational results, this does not 
happen for the age of the universe whose minimum value here results equal to 2/H . 
Let us finally treat the case which corresponds to the choice = 1 in the bulk. The relevant 
Einstein equations are 



alnl 



2 do^o 2 do 



a nl 



2 A, 



a nl 



As 
6 



1 /A 



sin 




+ 



1 /A 



A V 6 



cos^ 



AVIV 



6 a 2 . o 1 
sm 2 - 



3k 



3d 2 

2 2 



a 5 

2A 5 



A 
/ 




6 



As 
6 



sm z - 



A 



/A. ( 
6 



+ 



cos 2 — 



1 

A 




3k 



6 a 5 



sin 



A 




The solutions to the system (50) are 

/ 



c\ exp 



of 



^ 3 (k + a 2 ) cos 2 — 



A. 



7 



sin 



A 



2-k/(k+a 2 ) 



cos 2 — 



-A, 



a 2 /(3(fc+a 2 )) 



A2 



n = 



y (A 4 + K? (pi + p 2 + p 3 )) - k 



Comparing eqs. (50) and (9) we have defined the 4-D cosmological constant as A 4 = 



(50a) 



(50b) 



(51) 



(52) 



2 A, 



and have considered the remaining terms as the superposition of three different fluids with 
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pressures and densities given by 



2 3A 4 
Pi = -^Pu Pi= , — r- (53) 

P2 = - P2, P2 = " , — (54) 

1 / — A4 



4«| cos 2 -^/ — t 
1 3/c 

P3 = -oP3, P3 = -— ~ 2 " , — (55) 

3 4a 2 2 . 2 1 /-A 4 



«4 sin — \ I — : — t 



A 4 4 A V 4 

Again we notice that in the actual case it is neither possible to recover the brane metric 
evaluating the bulk metric at y = nor in general to give the explicit dependence of the 
scale factor on the proper time r even assuming k — 0. As a qualitative estimate one can say 
that the scale factor 0%, starting from zero, either increases indefinitely or, after increasing, 
reaches again the zero depending on the values of the parameters appearing in eq. (51). 



6. Conclusions 



We have studied brane-world cosmologies where the Einstein-Hilbert action is modified 
by a 4-D scalar curvature from induced gravity on the brane. To investigate cosmological 
solutions in the brane-bulk system, we extended the static Randall-Sundrum line element 
to a simple dynamical case and from the knowledge of the metric so obtained we solved the 
4-D Einstein equations on the brane. The general features of the model are presented and 
the evolution of physically meaningfull quantities can be determined choosing suitably the 
model parameters. In our description the cosmological fluid appears as a mixture of perfect 
fluids obeying simple equations of state with constant barotropic factors. Of course our 
model has to be implemented to can describe the distinct periods of the universe evolution, 
a task which will require to consider more complicated equations of state for real fluids. If 
one adopts an approach of this kind, the simple model described here may serve as a basis 
for obtaining more detailed braneworld solutions and, therefore, a better comparison with 
the accumulated cosmological observations. 
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